Weak ferromagnetism and other instabilities of the two-dimensional t-t' Hubbard 

model at Van Hove fillings 
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We investigate magnetic and superconducting instabilities of the two-dimensional t-t' Hubbard 
model on a square lattice at Van Hove densities from weak to intermediate coupling by means of the 
Two Particle Self-Consistent approach (TPSC). We find that as the next- nearest-neighbor hopping 
\t'\ increases from zero, the leading instability is towards an incommensurate spin-density wave whose 
wave vector moves slowly away from (7r,7r). For intermediate values of \t'\, the leading instability 
is towards d^2„y2-wave superconductivity. For larger \t'\ > 0.33t, there are signs of a crossover 
to ferromagnetism at extremely low temperatures. The suppression of the crossover temperature is 
driven by Kanamori screening that strongly renormalizes the effective interaction and also causes the 
crossover temperature to depend only weakly on t' . Electronic self-energy effects for large \t'\ lead to 
considerable reduction of the zero-energy single-particle spectral weight beginning at temperatures 
as high as T < O.lt, an effect that may be detrimental to the existence of a ferromagnetic ground 
state at weak coupling. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.10.Lp 
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I. INTRODUCTION 

Historically, the single-band Hubbard model was sug- 
gested independently by Gutzwiller JJ, Hubbard and 
Kanamori Q to gain insight into the origin of metallic 
ferromagnetism. However, despite enormous efforts 
that were undertaken to find answers to this question, 
only a few reliable results have been obtained even for 
this simplest possible microscopic model. The Hubbard 
model also exhibits a variety of other competing phases, 
including antiferromagnetic and superconducting phases. 

The first exact results for ferromagnetism were ob- 
tained in the strong coupling limit, t/ ^ oo , by Na- 
gaoka and Thouless ^ who showed that the ground 
state of the Hubbard model with one hole or electron 
is ferromagnetic at an infinitely large Coulomb repul- 
sion. That result did not answer the question of sta- 
bility to a finite concentration of holes in the thermo- 
dynamic limit. Improved bounds for the Nagaoka state 
have recently been derived f?| for various lattices in two 
and three dimensions. Ferromagnetic ground states also 
occur if one of the several bands of the model is disper- 
sionless (so-called Lieb's ferrimagnetism HI and fiat-band 
ferromagnetism 9]). Mielke and Tasaki proved the local 
stability of ferromagnetic ground states in the Hubbard 
model with nearly flat 0| and partially filled jllj bands. 
Ref. ^3 contains a short review of these works as well as 
new results for Hubbard models without the singularities 
associated with flat bands. A review of results 13] ob- 
tained for the simple one-band Hubbard model in the last 
few years as well as the results of Mielke and Tasaki sug- 
gest that the important ingredients for ferromagnetism in 
that model are (a) an interaction strength that is in the 
intermediate to strong coupling regime and (b) a band 



that exhibits a strong asymmetry and a large density of 
states near the Fermi energy or near one of the band 
edges. Metallic ferromagnetism at weak coupling, usu- 
ally known as Stoner ferromagnetism, has in fact been 
ruled out a long time ago by Kanamori based on 
the argument that the renormalization of the interaction 
strength brought about by T— matrix effects (Kanamori 
screening) would never allow the Stoner criterion to be 
satisfied when the density of states at the Fermi level 
p (Ep) is non singular. Physically, the largest possible 
effective interaction, according to Kanamori, is equal to 
the kinetic energy cost for making the two-particle wave 
function vanish when the two particles are at the same 
site. That energy scales like the bandwidth p[Ep)~^ so 
that the Stoner criterion 1 — Up [Ep] = cannot be ful- 
filled. Quantum Monte Carlo calculations confirm the 
quantitative nature of Kanamori's T— matrix result |l4j . 

If there is Stoner-type ferromagnetism in weak to in- 
termediate coupling, it is thus clear that, as in the mod- 
erate to strong-coupling case, one needs at least a sin- 
gular density of states to overcome Kanamori screening. 
An example of a model with singular density of states 
at the Fermi energy as well as band asymmetry is the 
two-dimensional (2D) Hubbard model with both near- 
est neighbor, t, and next-nearest-neighbor, t' , hoppings. 
When the Fermi energy is close to the Van Hove singu- 
larity the corresponding filling is usually referred to as 
a "Van Hove filling". At that filling, the Fermi surface 
passes through the saddle points of the single-particle 
dispersion. There are, however, other phases competing 
with ferromagnetism. At weak to moderate values of the 
on-site Coulomb repulsion U, for small t'/t and close to 
half-filling, the 2D t — t' Hubbard model shows an anti- 
ferromagnetic instability. That instability due to nesting 
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is however destroyed |l3 for a sufficiently large ratio t' /t 
at weak interactions in two and three dimensions, thus 
leaving room for other instabilities, including d-wave su- 
perconductivity and metallic ferromagnetism. 

The questions which we address in this paper are 
thus the following. Can the asymmetry of the band 
and the large density of states near the Fermi energy 
overcome the Kanamori argument and lead to ferromag- 
netism in the 2D Hubbard model? What are the com- 
peting phases? Most results on this problem (particu- 
larly for a square lattice) fall into three different classes, 
(a) Momentum-cutoff renormalization group (RG) meth- 
ods ITfil IT^ , and Quantum Monte Carlo calculations fl^ 
suggest that there is no evidence for ferromagnetism. 
But the problem, in particular with numerical meth- 
ods, is that only very small system sizes can be used 
in a regime where the size dependence is important. In 
addition, momentum-cutoff RG does not allow the con- 
tribution of ferromagnetic fluctuations these 
results should not be considered conclusive, (b) The 
second class of results is based on Wegner's flow equa- 
tions. They show j2fl] a- tendency towards weak ferro- 
magnetism with s*-wave character (the order parame- 
ter changes sign close to the Fermi energy). According 
to the flow equations calculations this phase competes 
with other instabilities in the particle-hole channel, in 
particular with the Pomeranchuk instability. The dif- 
ficulty of those weak-coupling calculations is that the 
s*-magnetic phase occurs at stronger coupling than the 
regime of validity of the second order analysis in U of 
the flow equations, (c) The third class suggests clear 
evidence for ferromagnetic ground states. These works 
include a projector Quantum Monte Carlo calculation 
with 20 X 20 sites and the T-matrix technique [23|, a 
generalized random phase approximation (RPA) includ- 
ing particle-particle scattering 22\ and exact diagonal- 
izations |23j |. Similar tendencies have been found by 
the authors of Refs. |3i US within the renormalization 
group and parquet approaches for the so-called two-patch 
model. Honcrkamp and Salmhofer recently studied [l^ 
the stability of this ferromagnetic region at finite temper- 
atures by means of a Temperature Cutoff Renormaliza- 
tion Group (TCRG) technique analogous to that used 
earlier for one-dimensional systems |2g. For U = 3, 
they have found that the ferromagnetic instability is the 
leading one for \t'\ > 0.33 |t| at Van Hove fillings with 
the critical temperature strongly dependent on the value 
of t' . When the electron concentration deviates slightly 
away from the Van Hove filling, the tendency towards 
ferromagnetism is cut off at low temperatures and a 
triplet p-wave superconducting phase dominates. The U- 
dependence of these ferromagnetic and superconducting 
phases in the ground state has been studied in Ref. [23 
by means of the same TCRG at weak coupling. 

In the present paper we study ferromagnetism and 
competing phases in the t — t' Hubbard model at weak to 
intermediate coupling by means of the two-particle self- 
consistent (TPSC) approach |2^. Antiferromagnetism 



and (i2;2_j^2-wave superconductivity are the competing 
instabilities. The TPSC approach is non-perturbative 
and applies up to intermediate coupling. It enforces 
the Pauli principle, conservation laws and includes the 
Kanamori screening effect. Comparisons with Quantum 
Monte Carlo calculations have shown that TPSC is the 
analytical approach that gives the most accurate results 
for the spin structure factor [29j |. the spin susceptibil- 
ity :28:| and the dj;2_j^2-wave susceptibility .'iO in two di- 
mensions. Throughout the paper we consider the 2D t—t' 
Hubbard model at Van Hove fillings from weak to mod- 
erate couplings. We determine the regions of the T — t' 
plane where the uniform paramagnetic phase becomes 
unstable to various types of fluctuations. We also esti- 
mate the electronic self-energy effects for large t' where 
ferromagnetic effects are present. The next section re- 
calls the methodology. We then present the results and 
conclude. 



II. TWO-PARTICLE SELF-CONSISTENT 
APPROACH 

We consider the t — t' Hubbard model on a square 
lattice with nearest (t) and next-nearest {t') neighbor 
hoppings 

+ 11^^71,^71,1, (1) 
i 

where c\^{cia) is the creation (annihilation) operator for 
the electrons with spin projection a € {j, |}, U is the lo- 
cal Coulomb repulsion for two electrons of opposite spins 
on the same site, and ni„ — cl^Cicr is the occupation num- 
ber. The bare single particle dispersion has the form, in 
units where lattice spacing is unity, 

ek = — 2t(cos kx + cos ky) — At' cos kx cos ky. (2) 

This spectrum leads to a Van Hove singularity in the den- 
sity of states coming from saddle points of the dispersion 
relation that are located at k = (0, ±7r) and (±7r, 0). The 
corresponding energy is Evh — 4i'. In this paper we al- 
ways consider the case where the non-interacting chem- 
ical potential is At', so that the non-interacting Fermi 
surface crosses the saddle points and the non-interacting 
density of states diverges logarithmically at the Fermi 
energy. The filling corresponding to this choice of chem- 
ical potential is a "Van Hove filling". For t' — and 
half-filling the Fermi surface is perfectly nested, namely 
£k+Q — ~£k, with Q — (tTjTt), which leads to an anti- 
ferromagnetic instability for J7 > 0. The perfect nesting 
is removed for t' /t ^ 0. We work in units where Bolz- 
mann's constant ks and nearest-neighbor hopping t are 
all unity. 
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The TPSC approach [23 can be summarized as fol- 
lows [S^- We use the functional method of Schwinger- 
Martin-Kadanoff-Baym with source field to first gener- 
ate exact equations for the self-energy E and response 
(four-point) functions for spin and charge excitations 
(spin-spin and density-density correlation functions). In 
such a scheme, spin and charge dynamical susceptibilities 
can be obtained from the functional derivatives of the 
source dependent propagator G with respect to 4>. Our 
non-perturbative approach then consists in two steps. 

At the first level of approximation, we use the follow- 
ing two-particle self-consistent scheme to determine the 
two-particle quantities: We apply a Hartree-Fock type 
factorization of the four-point response function that de- 
fines the product SG but we also impose the important 
additional constraint that the factorization is exact when 
all space-time coordinates of the four-point function coin- 
cide. From the corresponding self-energy, we obtain the 
local momentum- and frequency-independent irreducible 
particle-hole vertex appropriate for the spin response 
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The renormalization of this vertex mainly comes 0, 
from Kanamori screening Q. The double occupancy 
entering this equation is then obtained self- 
consistcntly using the fluctuation-dissipation theorem 
and the Pauli principle. More specifically, the Pauli prin- 
ciple, (n^) — (tIct), implies that 

while the fluctuation-dissipation theorem leads to an 
equality between the equal-time equal-position correla- 
tion ((n| — n|)^) and the corresponding susceptibility, 
namely 



(("T-"i)') 



n-2{n^ni), (4) 



where, using the short-hand q = (q, 2i7rmT), the summa- 
tion is over all wave vectors and all Matsubara frequen- 
cies with T the temperature, n the electron filling, and N 
the number of lattice sites. The latter equation is a self- 
consistent equation for the double occupancy, or equiva- 
lently for Ugp in Eq. (^J, because the spin-susceptibility 
entering the above equation is 



xoiq) 



1 - ^UspXoiq)' 



(5) 



where xo{q) is the particle- hole irreducible susceptibility 
including the contribution from both spin components 



y (a)^—Y^ /(^k) - /(£k+q) 
^°^1> TV 2-^ 2nrmT - + ek+q ' 



(6) 



with /(e) the Fermi-Dirac distribution function. Eq. Q 
is also known as the local-moment sum rule. The Green 



functions at this first level of approximation, G^^\ con- 
tain a self-energy E^^^ that depends on double-occupancy 
but since this self-energy is momentum and frequency 
independent, it can be absorbed in the definition of the 
chemical potential. In the above then, G^^^ is the bare 
propagator and xo is the bare particle-hole susceptibility 
both evaluated with the non-interacting chemical poten- 
tial fiQ corresponding to the desired filling. The irre- 
ducible charge vertex Uch = |^ + |^ strictly speaking 
is not momentum and frequency- independent. Neverthe- 
less, assuming for simplicity that it is, it can be simply 
found by using the fluctuation-dissipation theorem for 
charge fluctuations and the Pauli principle. 



N ^ 



with 



xo(g) 

l + hUchXo{q)' 



(7) 



The spin and charge susceptibilities obtained from 
Eqs. © and Q satisfy conservation laws [IE IIS- This 
approach, that satisfies the Pauli principle by construc- 
tion, also satisfies the Mermin- Wagner theorem: There is 
no finite-temperature phase transition breaking a contin- 
uous symmetry. Nevertheless there is a crossover temper- 
ature below which the magnetic correlation length grows 
exponentially 28] until it reaches infinity at zero tem- 
perature. Detailed comparisons of the charge and spin 
structure factors, spin susceptibility and double occu- 
pancy obtained with the TPSC scheme are in quantita- 
tive agreement with Quantum Monte Carlo simulations 
for both the nearest-neighbor [2^ and next-nearest- 
neighbor [s^ l Hubbard model in two dimensions. 

In loop expansions, response functions are computed 
at the one-loop level and self-energy effects appear only 
at the two loop level. Similarly, in our case the second 
step of the approach gives a better approximation for 
the self-energy. We start from exact expressions for the 
self-energy with the fully reducible vertex expanded in 
either the longitudinal or transverse channels. These ex- 
act expressions are easy to obtain within the functional 
derivative formalism. We insert in those expressions the 
TPSC results obtained at the first step, namely Usp and 

Uch, xip (9), Xchiq) and G(^)(fc-|-g) so that Green func- 
tions, susceptibilities and irreducible vertices entering the 
self-energy expression are all at the same level of approx- 
imation. Then considering both longitudinal and trans- 
verse channels, and imposing crossing symmetry of the 
fully reducible vertex in the two particle-hole channels, 
the final self-energy formula reads |^, 33] 

E(2)(fc) = Un, + ^^J2[^UspxilHQ) 

q 

+ UcnX^^^{q)\G^^\k + q). (8) 
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This self-energy © satisfies [2^ |^ the consis- 
tency condition between single- and two-particle prop- 
erties, Tr(S(2)G(i)) = 2U{n^ni). Internal consistency 
of the approach may be checked by verifying by how 
much Tr(E(2)G'(2)) differs from 2U{n^ni). The results 
for single-particle properties given by the self-energy for- 
mula JHJ are in quantitative agreement [28l [^S*, '34] with 
numerical simulations at weak to moderate couplings. At 
temperatures much lower than the crossover temperature 
where the correlation length increases exponentially, the 
consistency condition signals that the method becomes 
less accurate, although it does extrapolate in most cases 
to a physically reasonable zero temperature limit . In 
the present paper, we will not present results below the 
crossover temperature so we are always within the do- 
main of validity. It should be noted that the self-energy 
Eq. (O takes into account the fluctuations that are dom- 
inant already at the Hartree-Fock level, namely the an- 
tiferromagnetic ones. 

The above formalism can be extended [30[ to com- 
pute pairing correlations. Physically, the dx2-y2-wsLve 
susceptibility shows up after antiferromagnetic fluctua- 
tions have built up since it is the latter that give some 
non-trivial momentum dependence to the vertices. Mo- 
mentum dependence of the vertices is absent in the bare 



Hamiltonian and also at the first level of TPSC. It ap- 
pears from the momentum dependence of the self-energy 
at the second level of approximation. In other words, our 
formalism physically reflects old ideas about pairing by 
antiferromagnetic spin waves [s^ . What it contains that 
is absent in other formalisms is the possibility of sup- 
pression of superconductivity by pseudogap effects also 
induced by antiferromagnetic fluctuations 30]. 

The mathematical procedure to obtain the dx2_y2- 
wave pairing susceptibility is as follows. Basically, the 
above steps are repeated in the presence of an in- 
finitesimal external pairing field that is eventually set 
to zero at the end of the calculation. This allows 
us to obtain the particle-particle irreducible vertex in 
Nambu space from the functional derivative of the off- 
diagonal S*^^^ with respect to the off-diagonal Green 
function. The d-wave susceptibility is defined by Xd — 
JI^ dr (Tt-A (t) A^"^ with the d^2_y2-wsive order param- 
eter = J2lJ2J9h)4^4+Jl' t^^c sum over 7 being 
over nearest- neighbors, with 5(7) — ±1/2 depending 
on whether 7 is a neighbor on the x or on the y axis. 
(3 = 1/T, Tr is the time-ordering operator, and r is 
imaginary time. The final expression for the dj.2_y2-wave 
susceptibility is 



Xd 



(q = 0, zg™ = 0) = ^Y1 ^ {-k) cf^ (k)) - 
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with gd(k.) = {cos k,j; — cos ky) the form factor appropri- 
ate for d-wave symmetry. The above expression contains 
only the first two-terms of the infinite series correspond- 
ing to the Bethe-Salpeter equation. It should be noted 
that the appearance of G^^^ on the right-hand side of the 
equation for the susceptibility Eq. ((UJ allows pseudogap 
effects to suppress superconductivity 30] . This effect is 
absent in conventional treatments of pairing induced by 
antiferromagnons . 

Since the crossover to the ferromagnetic ground state 
found in our work appears at very low temperatures 
(T < 1/200), a large lattice is required in order to avoid 
finite-size effects at those temperatures. In the case of fer- 
romagnetism, sensitivity of the results to the lattice size 
at low T can be avoided by making sure that the lattice is 
large enough at any given temperature to reproduce the 
weak In T behavior of the bare particle- hole susceptibility 
Xo (q =0, iqm = 0). That singularity reflects the singular 
density of states at the Van Hove filling. We found that a 
N = 2048 X 2048 lattice suffices to compute xo entering 
the TPSC phase diagram. The sum over q in Eq. Q 



can be performed on a coarser mesh without loss of pre- 
cision. To speed up the calculations and to overcome 
increasing memory requirements, especially at low tem- 
peratures, we use the renormalization group acceleration 
scheme .36] . Interpolation is used to obtain quantities at 
temperatures that fall between those directly obtained 
with the renormalization group acceleration scheme. 



III. WEAK FERROMAGNETISM AND OTHER 
INSTABILITIES 

Without loss of generality, we can take t > and t' < 0. 
In that case, the Van Hove filling is always at n < 1. The 
Van Hove fillings n > 1 occur only when t and i' have 
the same sign, but this case can be mapped back to the 
situation n < 1 using the particle-hole transformation 
— > {—lydia and Cia (~1)''^L where the phase 
factor takes the value -1-1 on one of the two sublattices 
of the bipartite lattice and —1 on the other sublattice. 
The sign of t and f can be changed simultaneously with 
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the particle-hole transformation defined by c]^ di^ 
and Cia d\^. Whenever a particle- hole transformation 
is performed, the occupation number changes from n to 
2 — n. The Van Hove filling vanishes at \t'\ = 0.5 jt| so 
we restrict ourselves to \t'\ < 0.5 For larger \t'\ there 
is a change in Fermi surface topology. 

We begin with the Random Phase Approximation 
(RPA) phase diagram in the T ~ t' plane, then move to 
the TPSC crossover diagram and conclude with a short 
section on effects that can be detrimental to ferromag- 
netism. 



A. RPA phase diagram 

Within RPA or mean-field, the transition temperature 
Tc may be found from 
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2-C/xo(q,0) = 0, 



(10) 



where XoCQjO) is the zero- frequency limit of the non- 
interacting particle-hole susceptibility given by Eq. ©. 
In the case of ferromagnetism q = (0,0), while q = 
Q = (tt, tt) in the case of commensurate antifcrromag- 
netism. The temperature at which the uniform para- 
magnetic phase becomes unstable to fluctuations at the 
antiferromagnetic or at the ferromagnetic wave vector is 
plotted in Fig. One should keep in mind that, in all 
cases, we are speaking of spin-density waves, namely the 
local moment is in general smaller than the full moment. 
Furthermore, for \t'\ different from zero, the real wave 
vector where the instability occurs is incommensurate. 
The question of incommensurability is considered in more 
details in the TPSC section. Note that in contrast to 
the case U = 3, the ferromagnetic critical temperature 
for U = 6 does not increase with t' , it even decreases 
slightly. We do not explore the stability of the various 
phases that could occur in mean-field theory below the 
indicated transition lines. 

In both RPA and TPSC, the wave vector where the 
instability first develops is related to the q-dependence 
of xo- In TPSC, it is not only the maximum value of 
Xo (QiO) that determines the crossover temperature, but 
the whole q-dependence of xo that comes in the sum rule 
Eq. 0) for Usp. From the plot of xo as a function of wave 
vector at T = 0.01 in Fig. El one can see that at t' = 
the antiferromagnetic wave- vector Q leads to the largest 
value of Xo- With increasing \ f\ the maximum of xo is 
at an incommensurate wave vector Q5 — {t: — S, tt) close 
to (tt, tt), while for large \t'\ > 0.32 the maximum moves 
clearly to (0,0). For intermediate negative values of the 
next-nearest-neighbor hopping \t'\ 0.3 the magnitudes 
of the susceptibility at (0, 0) and at (tt, tt) are comparable 
so the change in the relative magnitude as a function of 
temperature is important. 

The main deficiencies of RPA are (a) finite temper- 
ature phase transitions in two dimensions that contra- 
dict the Mermin- Wagner theorem, (b) an overestimation 
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FIG. 1: The RPA critical temperature Tc as a function of the 
Van Hove fillings indicated on the upper horizontal scale and 
the corresponding value of next-nearest-neighbor hopping t' 
on the lower horizontal scale. The critical temperature Tc is 
determined from Eq. 1)10^ . AFM stands for the region where 
the uniform paramagnetic phase becomes unstable to fluctu- 
ations at (tt, tt) while FM is the region where the instability 
is at (0, 0). Vertical lines denote the boundary between AFM 
and ferromagnetic phases. 



of the effect of U on Tc because of the neglect of the 
renormalization of U brought about by quantum fiuctua- 
tions (Kanamori screening). One can see from Fig.^that 
the RPA critical temperature is quite a bit larger than 
the crossover lines predicted by the TCRG (see Fig. 1 of 
Ref. ^3)- The TPSC remedies these deficiencies. 



B. TPSC crossover diagram 

We begin by considering the effective interaction Usp 
that plays a crucial role in TPSC. In Fig. O we plot 
Usp as a function of t' as obtained from Eqs. l|2Jl, Q 
and (O. One can see that Kanamori screening strongly 
renormalizes the effective interactions. This weakly tem- 
perature dependent renormalization effect is stronger for 
large |t'| in comparison with small \t'\. To explain this 
behavior we consider the sum rule that determines Usp, 
Eq. 10}. The main contribution to the sum on the 
left-hand side of this equation comes from the small 
denominator caused, for large \t'\ by xo(O,0), and for 
small \t'\ by Xo(QjO). As the coefficient before the log- 
arithm scales as [a/I — 4(t'/i)2]^i for xo(0,0), and as 

In {1 + ^l-4:{t'/ty)/{2t'/t) for xo(Q,0), it turns out 

that xo(O,0) increases rapidly for \t'\ near 0.5. This 
means that Usp has to decrease at large \t'\ to satisfy the 
sum rule Q where, in addition, the quantity n — 2(n|n|) 
on the right-hand side is a decreasing function of density 
(and hence of \t'\). 
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FIG. 2: The non-interacting particle-hole susceptibility xo at 
zero frequency as a function of wave vector q along a path in 
the Brillouin zone is drawn for various values of next-nearest- 
neighbor hopping t' a,t T = 0.01. The filling is obtained by 
placing the chemical potential at the energy of the Van Hove 
singularity for the given t' . 
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FIG. 3: Irreducible spin vertex Usp as a function of next- 
nearest-neighbor hopping t' (or corresponding Van Hove fill- 
ings on the upper horizontal scale) at T = 0.125. Horizontal 
lines at [/ = 3, 6 denote the bare Hubbard repulsion U. 



To find the crossover lines, we consider the zero- 
frequency limit of the spin susceptibility given by Eq. 
and the dx^-y^-wave pairing susceptibility given by 
Eq. © above. The crossover temperature Tx for the 
magnetic instabilities is chosen as the temperature vifhere 
the enhancement factor Xspi<l, 0)/xo(q, 0) is equal to 500. 
We have checked that this corresponds to a magnetic cor- 
relation length that fluctuates around 25 lattice spacings 
for \t'\ between \t'\ — and \t'\ — 0.3. The crossover 
temperature Tx is not very sensitive to the choice of cri- 
terion because near and below the crossover region the 
enhancement factor grows very rapidly (exponentially). 

For pairing, we proceed as follows. Eq. ((SJ contains 
only the first two terms of the infinite Bethe-Salpeter 
series. The first term (direct term) describes the prop- 
agation of dressed electrons that do not interact with 
each other while the second term contains one spin fluc- 
tuation (and charge fluctuation) exchange. This comes 
about in our formalism because S*-^) is a functional of 
G'"'^' . We would have obtained an infinite number of spin 
and charge fluctuation exchanges, in the usual Bcthe- 
Salpeter way, if we could have written S*^^^ as a functional 
of G^^). This is not possible within TPSC. We have only 
the first two terms of the full series. The superconduct- 
ing transition temperature in two dimensions is of the 
Kosterlitz-Thouless type and is expected to occur some- 
what below the temperature determined from the Bethe- 
Salpeter equation (Thouless criterion). We thus use, as a 
rough estimate for the transition temperature for d-wave 
superconductivity, the temperature where the contribu- 
tion of the vertex part (exchange of one spin and charge 
fluctuation) becomes equal to that of the direct part (first 
term) of the d-wave pairing susceptibility 30]. In other 
words, we look for the equality of the sign and the mag- 
nitude of the two terms appearing in Eq. @ . This choice 
is motivated by the statement that 1 -I- a; -f . . . resumnied 
to 1/ (1 — a;) diverges when a; = 1. 

The TPSC phase diagram shows three distinct regions 
illustrated for [/ = 3 and for U = 6 in Fig. 21 (a) for 
t' — 0, the leading instability is at the antiferromag- 
netic wave vector and for small non- vanishing \f\ it is 
at an incommensurate wave vector close to (7r,7r). We 
will loosely refer to that region as the region where anti- 
ferromagnetism dominates, (b) For intermediate values 
of the next-nearest-neighbor hopping, c?3.2_y2-wave super- 
conductivity dominates, (c) At large negative \t'\ > 0.33 
a crossover to a magnetic instability at the ferromag- 
netic wave vector occurs. Let us consider these different 
regions in turn. 

Near t' = 0, Tx is relatively high and the suscepti- 
bility near the antiferromagnetic wave vector grows most 
rapidly. When we increase \t'\, the crossover temperature 
decreases because of reduced nesting of the Fermi surface. 
In TPSC the wave vector of the instability is incommen- 
surate for any finite value of the next-nearest-neighbor 
hopping \t'\ as can be concluded from the structure of 
Eq. jSJ and from the fact that the non-interacting suscep- 
tibility with momenta = (ir — S, tt) is the largest when 
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asymptotic behavior of the Lindhard function near q 
and q = Q is, respectively, 21l 



FIG. 4: The TPSC phase diagram as a function of next- 
nearest-neighbor hopping t' (lower horizontal axis). The cor- 
responding Van Hove filling is indicated on the upper hori- 
zontal axis. Crossover lines for magnetic instabilities near the 
antiferromagnetic and ferromagnetic wave vectors are repre- 
sented by filled symbols while open symbols indicate insta- 
bility towards d^2_y2-wa,ve superconducting. The solid and 
dashed lines below the empty symbols show, respectively for 
U = 3 and U = 6, where the antiferromagnetic crossover 
temperature would have been in the absence of the super- 
conducting instability. The largest system size used for this 
calculation is 2048 x 2048. 



i' 7^ 0. The incommensurate wave vectors are plotted in 
Fig. [SI as a function of t' . One can see that the degree of 
incommensurability is strongly temperature-dependent, 
and that it increases with increasing temperature. 

In the second region of the TPSC phase diagram 
(ij.2_^2-wave superconductivity is the leading instability. 



Xo(O,0) ^ ln(l/max(Ai,r))/v/l-i?2 
Xo(Q,0) ^ ln(l/ max(A., T)) In [(1 + ^1 - R^)/R 



with R = 2t'/t so that, looking at the equality of 
the coefficients of the logarithms, one finds that the 
change from antiferromagnetic to ferromagnetic behav- 
ior occurs atjt'l = 0.27 instead of \t'\ ~ 0.33 as found 
above p^. l38l| . To understand the difference between 
these two results, we need to look at subdominant cor- 
rections. For example, a numerical fit reveals that 
Xo(Q,0) ~ 0.52 + 0.241ogio(l/T). This means that for 
the leading term with a logarithmic structure to be, say, 
about ten times larger than the subdominant term, the 
temperature should be as low as 10~^°. The correspond- 
ing U (or Usp) that satisfies 1 = [/(or C/sp)xo(Q, 0)/2 at 
this temperature is very small, namely OAt. Therefore, 
unless U is very small, the next leading term plays an 
important role and a straightforward application of the 
asymptotic form (taking only the leading term) is not 
justified. For U — 6 and U — 3, for example, TPSC 
shows that near the antiferromagnetic to ferromagnetic 
boundary, the crossover temperature is of order 10^^ and 
lO^'^ respectively. For this temperature, the sub-leading 
term 0.52 is comparable to the logarithmic contribution. 

The TPSC phase diagram is in qualitative agreement 
with the TCRG phase diagram In addition, the crit- 
ical values t'^ for the stability of superconductivity and 
ferromagnetism are the same in both approaches. But in 
contrast with the TCRG, ferromagnetism in TPSC oc- 
curs at very low temperatures, and increasing |t'| does 
not lead to a dramatic increase in crossover tempera- 
ture. One can see from Fig. 0] that the critical values of 
t' for the stability of ferromagnetism are unchanged for 
different U, whereas the critical \t'J for the stability of 



In this regime the transition temperature to ^2. 2 -wave 42_j^2-wave superconductivity decreases with increasing 



superconductivity appears higher than the temperature 
at which the antiferromagnetic correlation length be- 
comes larger than about 25. The latter crossover lines 
are denoted by the solid {U = 3) and by the dashed 
lines {U = 6) in Fig. 0| Note that da;2_j,2-wave super- 
conductivity is here induced by incommensurate antifer- 
romagnetic fluctuations. While high-temperature super- 
conductors are not generally close to Van-Hove singu- 
larities, incommensurate dynamic spin fluctuations are 
concomitant with d.j.2_y2 superconductivity in these com- 
pounds "37|. 

Finally, the third regime occurs at \t'\ > 0.33 where the 
ferromagnetic susceptibility Xsp{0, 0) is the leading one at 
low temperatures. Ferromagnetism occurs because of the 
diverging density of states at the Van Hove singularity. 

Note that for U infinitesimally small the phase bound- 
aries happen close to zero temperature. Disregarding su- 
perconductivity for the moment, let us consider where 
the phase boundary between antiferromagnetism and fer- 
romagnetism would be at small U. In that case, the 



coupling strength U. 

The fact that the crossover temperature towards ferro- 
magnetism depends even more weakly on t' in TPSC than 
in RPA can be explained by the following simple argu- 
ment. Taking into account Kanamori's improvement [Sj 
of the naive Stoner criterion for ferromagnetism, we ex- 
pect that the crossover temperature Tx can be roughly 
approximated by 



Tx ^ To exp 



1 



(11) 



where Tq is a constant, p{Ef) — xo(0,0)/2 and J/eff 
is the renormalized effective interaction [Ugp in the case 
of TPSC). We have already explained in the context of 
Fig. O that the increase with \t'\ of the weight of the log- 
arithmic singularity in the density of states at the Fermi 
level leads to a decrease of Usp , so the crossover temper- 
ature is almost constant in TPSC. 

A distinctive feature of the TPSC phase diagram is 
that the crossover to ferromagnetism generally occurs at 
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FIG. 5: Incommensurate wave vector — {n — S, n) where 
the maximum of the non-interacting susceptibility is located 
as a function of next-nearest-neighbor hopping t' at Van Hove 
fillings. Different lines correspond to different temperatures. 
Given t' and a crossover temperature in the TPSC phase di- 
agram, one can use the present figure to find out the incom- 
mensurate wave- vector at which the instability first occurs. 



much lower temperature than the crossover to antifer- 
romagnetism. This partially comes from the pecuhar- 
ity of the temperature dependence of the zero-frequency 
limit of the non-interacting particle-hole susceptibility. 
To demonstrate this, let us use, as an estimate for 
the crossover temperatures in TPSC, the RPA criterion 
Eq. H1U|) with U replaced by Usp and let us look for values 
of the temperature when the left-hand side of that equa- 
tion becomes small (it will vanish only at zero tempera- 
ture). At small \t'\ the leading non- interacting staggered 
susceptibility Xo(Q7 0) behaves like (InT)^ with temper- 
ature, while for |<'| > 0.33 the leading non-interacting 
uniform susceptibility xo{0,0) scales as |lnT|. We find 
that these susceptibilities have comparable size for tem- 
peratures r > 1, while the divergences of xo(Q,0) and 
Xo(0, 0) begin respectively at T < 1 and T <C 1. There- 
fore, since the Stoner criterion Eq. (|10|l is satisfied in 
RPA with bare t/ = 3,6 at temperatures T > 1, RPA 
shows the same temperature scale for ferromagnetism 
and antiferromagnetism. But in TPSC the strong renor- 
malization of the interaction strength Usp < U means 
that the crossover occurs for larger values of Xo(QiO) 
and Xo(0, 0), in a regime where they already have differ- 
ent scales since Xo(Qi 0) for small \t'\ starts to grow log- 
arithmically at much higher temperature than xo(0, 0) 
for large \t'\. Thus, the crossover to antiferromagnetism 
in TPSC occurs at much higher temperatures than the 
crossover to ferromagnetism. 

Another interesting feature of the TPSC phase dia- 
gram at t/ = 3 is that the crossover temperature for anti- 
ferromagnetism is of the same order of magnitude as that 
of the TCRG result of Ref . , whereas the crossover to 
ferromagnetism is at much lower temperature than that 



observed in the TCRG calculations. The naive explana- 
tion is as follows. Let us assume that the approximate 
mean-field like expression Eq. (|11|) for Tx has meaning 
both within TPSC and within TCRG except that Uefj 
has a different value in both approaches. Simple algebra 
then shows that the relation between the crossover tem- 
peratures for TPSC and TCRG in the ferromagnetically 
fluctuating regime is 



ttpsc 



T^TPSC I 



1-1/a 



V^FM 



with a 



r/TCRG /TT 



characterizing the different renor- 
malizations of U in both approaches. When a — 1, 
both crossover temperatures are equal. For a > 1 the 
TCRG value for Tx is larger than for TPSC while the 
reverse is true when a < 1. Using the numerical re- 
sult [33 for the TCRG effective interaction at C/ = 3 and 
\t'\ ~ 0.45 we have a = 1.4 - 1.8. Then, replacing Tq 
by the bandwidth 8t and taking Tj^^'^ = 3.4 x 10"^ 
corresponding to \t'\ > 0.42 we obtain the estimate 
^ 10-30. This agrees with the crossover 
temperatures extracted from the TPSC (Fig. ^ and the 
TCRG phase diagrams (Fig. 1 of Ref. Il9|). Similarly in 
the antiferromagnetically fluctuating regime near \t'\ — 0, 
we use the improved mean-field estimate for Tx 



Tx ^Toexp{~y/8t/U,« 



to extract the following relation between the crossover 
temperatures 



T^TCRG / 
^ AFM _ ( 
T^TPSC 
^AFM 



\ T^TPSC 

V-'afm 



Using the value of from the TPSC and the TCRG 
effective interaction 39] at U — 3 and |t'| ~ 0.1 we have 
a = 1.0 - 1.4. This leads to T^^^^ /T^^^^ « 1 - 2.5 for 
T^p^p ~ 4 X 10"^ at \t'\ ~ 0.1, which is in good agree- 
ment with the data extracted from the phase diagrams. 

As mentioned at the beginning of this subsection, the 
crossover temperatures Tx for the magnetic instabilities 
in TPSC have been chosen such that the enhancement 
factor is equal to 500. The enhancement factor scales 
like the square of the correlation length . For such large 

the value of Tx is rather insensitive to the choice 500 
since the correlation length grows exponentially. Our cri- 
terion for Tx leads to a good estimate of the real phase 
transition temperature with ^ = 00 when a very small 
coupling term is added in the third spatial direction. 
The dependence of Tc on coupling in the third dimen- 
sion has been studied, within TPSC, in Ref. EOl. The 
latter reference also contains expressions for the relation 
between the enhancement factor and On the other 
hand, Tx depends more strongly on the precise criterion 
if we choose a moderate value of the enhancement factor. 
In particular, the TPSC value of Tx in the antiferromag- 
netic fluctuation region increases by a factor two to five 
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if we choose 10 for the enhancement factor, close to the 
value ^3 chosen in Ref. p^]. In this case, Tx agrees es- 
sentially perfectly with the value obtained in the TCRG 
phase diagram. 

Note however that our estimate for the superconduct- 
ing transition temperature is smaller than that obtained 
with the TCRG of Ref. P). Because in TPSC the pair- 
ing fluctuations do not feed back in the antiferromag- 
netic fluctuations, this result suggests that the feedback, 
usually included in TCRG, enhances superconductivity 
in this region of the phase diagram. Such a positive 
feedback effect was also found in the RG calculations of 
Refs. ITtI fi^ . On the other hand, the RG approach of 
Ref. [2J| suggests instead that antiferromagnetism and 
superconductivity oppose each other. Some particle- 
particle diagrams were however neglected in the latter 
approach. In TPSC, antiferromagnetic fluctuations help 
d2,2_j^2-wave superconductivity as long as they are not 
strong enough to create a pseudogap, in which case they 
are detrimental to superconductivity |30| . 

The above-mentioned rcnormalization group calcula- 
tions were done in the one-loop approximation with- 
out self-energy effects. By contrast, in the RG work of 
Ref. 01 J self-energy effects showing up at two loops were 
included in the calculation for the t' = model. There, 
it was found that dressing the flow equations for AFM 
and superconducting response functions with the one- 
particle wave vector dependent weight factors Z results 
in a reduction of both AFM and superconducting cor- 
relations, the latter suppression being more pronounced. 
Within TPSC, the momentum- and frequency-dependent 
self-energy effects that appear in G(2) in the pairing sus- 
ceptibility Eq. lO do tend to decrease the tendency to 
pairing when AFM fluctuations become very strong at 
and near half-filling |3C|. in qualitative agreement with 
the RG result 01 • In particular, in the presence of an 
AFM-induced pseudogap, the tendency to superconduc- 
tivity is decreased compared to what it would be if we 
replaced G'^'^^ by G'^^' everywhere. (Such a replacement is 
not allowed within our formalism) . Because of the excel- 
lent agreement between TPSC at the first level of approx- 
imation and Quantum Monte Carlo calculations [281 12^ , 
momentum and frequency dependent self-energy effects 
are not expected to be very important for AFM fluctua- 
tions unless we are deep in the pseudogap regime. They 
have not been taken into account at this point. They 
might be more important in the case of ferromagnetism, 
which is already a very weak effect in TPSC. This is dis- 
cussed in the following subsection. 



C. Additional effects that may be detrimental to 
ferromagnetism 

The TCRG phase diagram is computed at the one- 
loop level. Self-energy effects occur at the two-loop level. 
Similarly, self-energy effects in TPSC are calculated at 
the second level of approximation. Since analytical con- 



T 




0.2 0.4 ^ 0.6 0.8 



FIG. 6: Temperature dependence of z'{T) defined by Eq. H12|l 
at the Van Hove filling corresponding to \t'\ = 0.4. 

tinuation of imaginary-time results is difficult at low tem- 
perature, we estimate the quasiparticle weight with the 
help of the quantity z'{T) defined in Refs. |2^0| by 

/(T).-2G(k.,;3/2)./ 1^4^54-. (12) 

J 27r cosh(pcL'/2) 

Physically, this quantity is an average of the single- 
particle spectral weight A{k.p,Lij) within T of the Fermi 
level {lo When quasiparticles exist, this is a good es- 

timate of the usual zero-temperature quasiparticle rcnor- 
malization factor z = (1 — dTi/duj)^^. However, in con- 
trast to the usual z, this quantity gives an estimate of the 
spectral weight A(k^,w) around the Fermi level, even if 
quasiparticles disappear and a pseudogap forms. 

Fig. shows the quasiparticle rcnormalization factor 
z' at a value \t'\ = 0.4 where ferromagnetic fluctuations 
dominate at very low temperatures. One observes a pro- 
gressive decrease of z' (T) with decreasing temperature. 
We checked that the single particle spectral function 
A(ki?,a;) begins to show a small pseudogap at the tem- 
perature where z' begins to decrease significantly. Since 
the ferromagnetic fluctuations are not yet strong enough 
at that temperature to create a pseudogap, this effect is 
completely driven by the singular density of states at the 
Van Hove flUing. In other words, second-order perturba- 
tion theory should suffice to observe the effect. The anal- 
ogous feature was previously found by one of the authors 
and his co-workers in a second-order perturbation 
study of the nearest-neighbor two-dimensional Hubbard 
model at half-filling. Self-energy effects near Van Hove 
points have also been discussed in Ref. [i^. The rather 
strong suppression of spectral weight at the Fermi wave 
vectors for temperatures larger than the crossover tem- 
perature found in the previous subsection would proba- 
bly reduce the true Tx or even completely eliminate the 
possibility of a ferromagnetic ground state if we could in- 
clude the feedback of this self-energy effect into the spin 
susceptibility. 
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The ferromagnetic fluctuation regime is also very sen- 
sitive to doping within TPSC. In fact, deviations of the 
filhng by 2 — 3% away from the Van Hove fiUing remove 
the crossover to the ferromagnetic regime. 

There is also an argument that suggests that a Stoner- 
type ferromagnetic ground state is unstable in the two- 
dimensional Hubbard model. Within RPA in the ferro- 
magnetic state |47j . the spin stiffness constant for spin 
waves in the ferromagnetic state is proportional to mi- 
nus the second derivative of the density of states at the 
Fermi level . Since the density of states as a function 
of energy (away from the Van Hove filling) has a positive 
curvature in two dimensions, that leads to a negative spin 
stiffness constant and thus to an instability. This argu- 
ment is based on the non- interacting density of states. 
The pseudogap effect mentioned in the previous para- 
graph changes the curvature of the density of states at 
the Fermi level and may stabilize the ferromagnetic state. 



IV. CONCLUSIONS 

As found within temperature-cutoff renormalization 
group (TCRG) [11 E^, TPSC suggests that ferromag- 
netism may appear in the phase diagram of the 2D t ~ t' 
Hubbard model at Van Hove fillings for weak to inter- 
mediate coupling. It is striking that the overall phase 
diagrams of TCRG and TPSC have some close similari- 
ties. As in TCRG, we find, Fig.01 that for small negative 
values of the next-nearest-neighbor hopping the leading 
instability is a spin-density wave with slightly incommen- 
surate antiferromagnetic wave vector (Fig. . We could 
find incommensurability at small \t'\ only for very large 
lattice sizes. The TCRG seems to indicate that very close 
to l^'l = 0, the wave vector remains pinned at (tt, tt) |23| 
but that could be due to the fact that coupling constants 
in TCRG represent a finite region in wave vector space 
and hence very small incommensurabilities cannot be re- 
solved. For intermediate values of we also find d^^-y^- 
wave superconductivity. The precise value of \t'\ for the 
onset of dx2_y2--wsLve superconductivity depends some- 
what on the criterion used for the crossover temperature. 
One clear difference with TCRG, however, is that the 
range of \t'\ where superconductivity appears increases 
with U whereas it decreases with U in TCRG [23|. At 
large \t'\ > 0.33t, a crossover to ferromagnetism occurs as 
a result of the diverging density of states. TPSC cannot 
tell us what happens below the crossover temperature, 
but that temperature is the relevant one in practice since 
any small coupling in the perpendicular direction would 
lead to a real phase transition. 

The critical value for ferromagnetism, \t'\ = 0.33t, co- 
incides with that found in TCRG JJ, .27]. This value of 
\f\ is smaller than that found within the T— matrix ap- 
proximation |2ll |. but that may be because of the cutoff 
to the Van Hove singularity imposed by the small sys- 
tem sizes used in that approach. The critical value for 
ferromagnetism, \t'\ = 0.33t, also differs from the value 



l^'l = 0.27t obtained in Ref. 01 in the limit of zero tem- 
perature. We have explained in Sec. IIII Ij| that for the 
crossover to occur sufficiently close to T = for the ar- 
guments of Ref. 24] to be correct, one needs values of 
U that are unrealistically small. At finite U (we studied 
[/ = 3 and U ^ 6), subdominant corrections to the log- 
arithms shift the critical \f/t\ = 0.27 found by Ref. |2J| 
to the value \t'/t\ = 0.33 found by us and TCRG. 

The differences between TCRG and other approaches, 
as well as their strengths and weaknesses, are well ex- 
plained in Refs. ulSjilJ- The smaller temperature scale 
for crossover to dj.2_y2-wave superconductivity in TPSC 
is a noteworthy difference between our approach and 
TCRG 01 . This may be due to the fact that our calcu- 
lations include self-energy effects which are absent in 
one-loop TCRG. But the most striking difference is the 
temperature scale for ferromagnetism that in our case re- 
mains extremely small away from the critical |<'| — 0.33i. 

We have shown that the low temperature scale for 
the crossover to ferromagnetic fluctuations comes from 
Kanamori screening that strongly renormalizes the effec- 
tive interaction (this effect is smaller in the antiferromag- 
netic regime). In TPSC this renormalization comes from 
the constraint that the spin response function with Usp 
should satisfy the local moment sum rule, Eq. This 
causes the crossover temperature to ferromagnetic fluctu- 
ations to depend weakly on t' and to remain small. As in 
the T— matrix approximation |2l|. Kanamori screening 
seems much stronger than what is obtained with TCRG. 
The latter approach perhaps does not include all the large 
wave vectors and large energies entering the screening of 
the effective interaction. 

Within TPSC then, the tendency to ferromagnetism 
seems very fragile. In addition, we checked that in 
TPSC ferromagnetism disappears for electron concen- 
trations that are only very slightly (2 — 3%) away from 
Van Hove fillin gs, in overall agreement with the results 
of the TCRG [la, So the question of the exis- 

tence of Stoner-type ferromagnetism at weak to inter- 
mediate coupling is not completely settled yet, despite 
the positive signs and the concordance of the most reli- 
able approaches. We have estimated the electronic self- 
energy effects for large \t'\ and found that the quasiparti- 
cle renormalization factor is reduced significantly at tem- 
peratures T < 0.1. As a result, the single-particle spec- 
tral function ACkp^uj) starts to show a small pseudogap 
which, at high temperature, is completely driven by the 
singular density of states, and not by the ferromagnetic 
fluctuations that appear only at very low temperature. 
This rather strong suppression of spectral weight at the 
Fermi wave vectors for T > Tx may further reduce Tx 
or even completely eliminate the crossover to a ferro- 
magnetic ground state. We have argued in Sec. IIII CI 
that other factors could be detrimental to a ferromag- 
netic ground state in two dimensions. In particular, as 
is the case with RG calculations [T^ l27| |. a consistent 
treatment of the electronic self-energy effects on the spin 
response function remains an open issue. 
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Another interesting problem for future investigations 
is the question of whether ferromagnetism could com- 
pete with the Pomeranchuk instability, i.e. a spon- 
taneous deformation of the Fermi surface reducing its 
symmetry from the tetragonal to the orthorhombic one. 
Temperature cutoff RG [23, |43| disagrees with a sugges- 
tion [l^ I20I that this is one of the possible leading 
instabilities of the 2D t — t' Hubbard model at Van Hove 
fillings. 

Note added in proof: B. Binz, D. Baeriswyl and 
B. DouQot [Ann. Phys. (Leipzig) 12, (2003); 
|cond-mat/0309645 have recently questioned the appli- 
cation of one-loop renormalization group to ferromag- 
netism, suggesting that the error produced by the one- 
loop approximation is of the same order as the term 
which produces the ferromagnetic instability. 
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